Quasi-equilibrium in glassy dynamics: a liquid theory approach 
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We introduce a quasi-equilibrium formalism in the theory of liquids in order to obtain a set 
of coarse grained long time dynamical equations for the two point density correlation functions. 
Our scheme allows to use typical approximations devised for equilibrium to study long time glassy 
dynamics. We study the Hypernetted Chain (HNC) approximation and a recent closure scheme by 
Szamel. In both cases we get dynamical equations that have the structure of the Mode-Coupling 
(MCT) equations in the long time region. The HNC approach, that was so far used to get equilibrium 
quantities is thus generalized to a fully consistent scheme where long-time dynamic quantities can 
also be computed. In the context of this approximation we get an asymptotic description of both 
equilibrium glassy dynamics at high temperature and of aging dynamics at low temperature. The 
Szamel approximation on the other hand is shown to lead to the exact Mode Coupling equation 
of Gotze for equilibrium dynamics. We clarify the way phase space is sampled according to MCT 
during dynamical relaxation. 



A longly debated issue in the theory of glassy relax- 
ation of low temperature liquids is the relation between 
dynamical slowing down and equilibrium free-energy 
landscape. In schematic Mean-Field models, where both 
equilibrium and dynamics can be studied exactly, dynam- 
ical arrest and equilibrium ergodicity breaking appear as 
different sides of the same coin [1]. This observation is 
at the basis of the "Random First Order Transition The- 
ory" of the glass transition. Unfortunately, beyond these 
schematic models, when one tries to include liquid theory, 
the techniques and approximations involved in the study 
of equilibrium and dynamics differ considerably even at 
the mean field level. At the dynamical level one has Mode 
Coupling Theory (MCT) which provides an approxima- 
tion for an exact evolution equation of the density-density 
correlator and requires the static structure function as 
input [2]. For equilibrium one uses the replica method 
in conjunction with approximations like the Hypernet- 
ted Chain (HNC) partial rcsummation of virial series [3] 
or the small cage expansion [4]. Both approaches dif- 
fer on specific predictions, but agree on the same gen- 
eral picture. Ergodicity breaking appears to be driven 
by critical fluctuations (dynamical heterogeneities) with 
the same universal characteristics and growth of correla- 
tions. In recent times there has been considerable effort 
to reconcile the approaches in the "/? regime" , where, af- 
ter a first fast relaxation, correlation functions remain 
close to a plateau value. On the one hand, Szamel [5] 
has derived the MCT equation for the non-ergodicity pa- 
rameter using a particular closure of the replicated YBG 
hierarchy for the density correlation functions. On the 
other, an important dynamical quantity, the celebrated 
exponent parameter of MCT, has been related to ampli- 
tude ratios between equilibrium correlation functions [6] 
and has been computed in the framework of replicated 



liquid theory [7, 8]. Unfortunately, general approxima- 
tion frameworks capable to describe in a unified way both 
equilibrium and dynamics are at present lacking. One of 
the difficulties lies in the necessity of including conserva- 
tion laws (mass, momentum and energy) in the dynam- 
ical description at short times. However this difficulty 
should disappear in a long time "a regime" where correla- 
tions decay below their plateau value. In that case quasi- 
stationarity should hold in the dynamical equations. 

In this paper we show that progress can be achieved us- 
ing a coarse grained description where the details of the 
short times are neglected and one concentrates on the 
slow part of the dynamics. Despite the fact that glasses 
are out of equilibrium, emerging properties of glassy re- 
laxation, as temperature-time superposition principle [2] 
or the appearance of effective temperatures in the aging 
regime [9], can be rationalized through the hypothesis 
that quasi-equilibrium principles drive the dynamics dur- 
ing the phase space exploration in the long time regime 
[10]. In a nutshell, the two steps glassy relaxation at 
low temperature can be described as a first rapid decay 
inside a metastable state, followed by a very slow relax- 
ation in a complex landscape of metastable states in the 
long time a region. The short time regime is sensitive to 
the details of the dynamical rules of evolution, while the 
long time one is expected to be independent from such 
details and it can be thought as a random walk in phase 
space where the states that are available at a given time 
are selected according to their Boltzmann weight. 

This insight has been formalized in a recent paper 
[11] through the introduction of an appropriate Markov 
chain construction [12, 13]. Consider a system with 
time dependent configurations specified by S T and sub- 
jected to a Hamiltonian H(S) and an overlap function 
Q(S,S') = J dxdyw(x — y)ps(x)PS'(y), that measures 
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the similarity between density profiles of couples of con- 
figurations. The function w doesn't affect what follows, 
provided it is short ranged. One can define the Boltz- 
mann pseudodynamics for the system as the following 
Markov Chain 

M(S T+1 \S T ) = 6 5(Q(S T+1 , S T ) - C(t + l,r))(l) 

which implements the previous idea, assuming that the 
states available at a given time r + 1 are the ones that 
lie at an overlap C(r + l,r) from the configuration at 
time r. For mean field spin glasses, Boltzmann pseu- 
dodynamics provides an exact coarse grained version of 
slow dynamics when an infinitely long chain is considered 
if the value of C(r + 1, r) is such that the constraints do 
not affect the free-energy of the system [11]. This result 
is remarkable at the fundamental and practical levels. 
At the fundamental level it provides a basic probabilistic 
description of how configuration space is explored during 
glassy relaxation. At the practical level, it opens the pos- 
sibility to study long time dynamics through equilibrium 
techniques. 

The scope of this letter is to study Boltzmann pseu- 
dodynamics in glassy liquids, and to discuss its relation 
with present dynamical theories, in particular, with the 
Mode Coupling Theory (MCT) of the glass transition and 
the theory of aging dynamics as is known from schematic 
spin models. As usual, we will need approximations to 
deal with the strongly interacting character of the sys- 
tems. We will concentrate in liquid theory approximation 
leading to closed equations for density-density correlation 
functions: 1) the well known Hypcrnctted Chain approx- 
imation 2) the Szamel closure scheme of the replicated 
YBG hierarchy [5]. Our main results are the following: 

• Within the HNC approximation we get a fully self- 
consistent dynamical equation for correlation and 
response functions that consistently describe the 
evolution of liquids at low temperature. 

• Starting from an equilibrium-initial condition our 
equations reduce to a single one for the correlation 
function. This equation displays a MC singular- 
ity at a previously identified dynamical transition 
Td- All universal properties of Mode Coupling The- 
ory are found within our approach. In particular, 
relaxation in the early and late j3 regime is char- 
acterized by power laws whose exponents a and 
b verify the MC relation L(l - a) 2 /T(l - 2a) = 
T(l + b) 2 /T(l + 2b) = A. The exponent parameter 
A coincides with the one computed in [7, 8] with 
completely different techniques. Differently from 
MCT which needs the external input of the static 
structure factor, we get here a complete set of equa- 
tions determining static and long time dynamical 
quantities. 



• Starting from a strongly off-equilibrium situation 
we find an aging regime. Modulo the identifica- 
tion of a soft mode of slow relaxation, the gross 
picture follows the main lines found in mean field 
spin glasses [14]. The system falls asymptotically 
in a non-equilibrium state where time translation 
invariance and fluctuation dissipation theorem are 
violated. We compute self-consistently the emerg- 
ing effective temperature. 

• As an alternative approximation, we use a scheme 
recently proposed by Szamel. We show that this 
approximation allows to derive the canonical Mode 
Coupling equation for the slow part of the momen- 
tum dependent correlator in the critical regime. We 
thus unveil that, as described by MCT, the re- 
laxation in the whole a regime consists in quasi- 
cquilibrium exploration of the configuration space 
and explain the insensitivity of this regime from the 
details of dynamical evolution [15, 16]. 

Boltzmann pseudodynamics and replicas. In order to 
study the pseudodynamic process one can start from the 
probability of a trajectory 

e -f3 H[S ] 

P(Sl,S L -i, So) = — ~zl$) — ( 2 ) 

L-l -/3H[St+i] 
' n Z[S r MT + l,T)] S iQ{S ^ S ^ - + 

where Z(p) = £s e-,3ff[S1 and Z[S, C] = 
J2 S > e~P H ^5(Q(S',S) - C) are normalization factors. 
The choice of /3 = /? corresponds to choosing an equilib- 
rium initial condition at temperature /3, while /3q = is 
useful to model a high temperature quench and an aging 
situation. The probability (2) can be used to compute 
averages of observables or the free-energy of the system 
as a function of time. To use it in practice, the replica 
method can be employed to deal with the partition 
functions Z[S T , C(r + l,r)] in the denominator through 
the substitution g , g c\t+i tJ] Z[S T , C(r+l,r)] n-r_1 . 
As usual, the numbers of replicas n T1 that can here 
depend on time, is kept integer for intermediate com- 
putations and send to zero at the end. Enforcing the 
constraint in (1) through a Lagrange multiplier v(t), 
and expliciting the partition functions as sums over 
configurations, the replicated system formally describes 
the equilibrium of an interacting mixture of different 
particle species (one for each time and replica index) and 
static approximation schemes can be applied. In [11] 
schematic long range spin glass models were analyzed. 
It was shown that, the resulting mean field equations 
can be identified with the dynamical equations in the a 
regime if the value of C(r + l,r) is itself chosen by the 
dynamics, and the corresponding Lagrange multiplier 
v{t) vanishes. 
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Here we want to see what happens if we move from the 
schematic model situation to the replicated liquid theory 
framework. 

Pseudodynamics in the HNC approximation. Let us 
introduce the basic fields that we want to study. Suppose 
that we start with a chain of systems where each system 
is composed by N particles interacting through a pair 
potential 4>{x). 

We can introduce the replicated one and two points 
density functions. Defining for notational transparency 
Greek indexes that code for both time and replica in- 
dexes, e.g. a = (r, a), 

N 

p a (x)=Y / S(x-4 a) ) (3) 

p a p(x, y) = (p a (x)pp(y)) - (p a (x))S a ,f}5(x - y) . (4) 

The starting point of our analysis are the replicated 
HNC equations 

ln[h a i3(x,y) + 1] + f34> al3 (x,y) = h a p{x,y) - c a p{x,y) 

(5) 

where 4> a p(x,y) = 4>aT;ba{x,y) = ^T^a,b<t>{x — V) + 
o~b.io~T,cr+iw(x — y)v(a) contains the inter particle poten- 
tial as well as the constraint, p 2 h a p{x, y) = p a p,(x, y) — p 2 
and c a /3(x,y) is the direct correlation function defined 
through the replicated Ornstcin-Zcrnike equation 

c a p(x, y) = h a p(x,y)-p^2 / dzh a ^(x,z)cjp(z,y), (6) 

7 J 

where we have supposed that the density (p a (x)) = p is 
constant. Notice that for v(t) = 0, these equations are 
formally identical to the usual HNC equations for repli- 
cated liquids and admit well known time independent 
solutions that describe a glass transition at well defined 
temperature in the static replica formalism [3]. Here 
we want to go beyond these solutions to study slow dy- 
namics. As explained in [11] this can be done seeking 
for continuous solutions in the infinite chain limit, where 
the time variables t = t/L etc. become continuous. The 
simplest paramctrization allowing us to achieve this goal, 
writing explicitly time and replica indexes, reads: 

h a p(x, y) = h(x, y; s, u) + 5 ab 6{s - u)Ah{x, y; s, u)+ 
+ Q(u — s)S a iAh(x, y; s, u)+ 
+ 6(s - u)SbiAh(x, y; s, u) 

(7) 

where here s and u are pseudotime indexes that tell the 
position along the Markov chain while a and b are replica 
indexes that index the specific replica at a fixed position 
along the chain. The function 9(x) is equal to one when 
x > and zero otherwise. An analogous expression can 
be written for the direct correlation function. 



In the limit of zero replicas n u 0, h(x, y; s, u) rep- 
resents the normalized density-density space and time 
dependent correlation function along the chain; the func- 
tion 

R h (x,y;s,u) = (3B(u- s) (8) 

cm 

for s > u is the response function that encodes how 
the density in x at pseudotime s varies when a small 
perturbing pressure is applied to the system in y at 
time u. As the effect of short time dynamics is con- 
densed in a single step of the chain, we can expect 
that h{x,y]s,u) is discontinuous at equal time with 
Ah(x, y; s, s) = h(x,y,s,s) — lim u _> s h{x, y; s, u) repre- 
senting the jump in correlation from the instantaneous 
to the plateau value. 

In what follows we will search for a space translational 
invariant solution of the equations. By using (8) and 
by taking the infinite chain limit we can derive from the 
Ornstein-Zernike equation in Fourier space the dynamical 
equations 

= W g [h,c] + ^ / dzh{q;s,z)R c (q;u 7 z) + 

p r (9) 

+ /? J dz Rh ( q > s > z "> c ( q > z > u "> 

dzR h {q;z,s)R c (q;u,z) (10) 

Ah(q; s, s) ~ Ac(q; s, s) + pAh{q; s, s)Ac(q; s, s) (11) 

where we have defined W q [h,c] = c(q;s,u) — 
h(q;s,u) + p[h(q; s, 0)c(g; u, 0) + h(q; s,u)Ac(q;u,u) + 

c(q; s, u)Ah(q; s, s)\ and W<j [Rh, R c ] = 
R c (q;u,s) - Rh(q;u,s) + pR h {q;u, s)Ac(q;u,u) + 
pR c (q;u, s)Ah(q; s, s). These equations are valid for 
u > s; the equivalent ones hold for s > u. Moreover 
from the HNC equations we obtain for s u 

ln[/i(x; s, u) + 1] = h(x; s, u) — c(x; s, u) (12) 

R c (x;s,u) = R h (x;s,u) fe ( x ' s ' u ) ( 13 ) 
n(x; s, u) + 1 

Notice the quantities Ah(q; s, s) and Ac(q; s, s) are actu- 
ally s-independent and represent the jump from the in- 
stantaneous to the plateau values of the correlation func- 
tions as effect of the short time dynamics that we neglect. 
The dynamical equations (9-12) provide a complete set of 
equations that can in principle be solved step by step in 
time. Moreover, they are reparamctrization invariant in 
the sense that if we know a solution h(q; s, «), Rh{q_] s, u), 
another solution can be obtained by putting h'(q, s, u) = 
h(qJ(s)J(u)) and R' h (q,s,u) = ^ R h (q; f (s) J (u) 
(with the analogous expressions for c and R c ) provided 
that f(x) is a monotonous increasing function such that 
/(0) - 0. 
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Equilibrium dynamics. At this point we are equipped 
to study the equilibrium dynamics, where one starts 
from an initial condition well thcrmalized at /3q = /3- 
For generic temperature, solutions with non-trivial time 
dependence, require non zero value of f(t), implying a 
change in the free-energy of the system. However, for 
T — > T£ a solution with vanishing v{t) emerges that sat- 
isfies the equilibrium properties of TTI and FDT. In this 
way the two equations (9) and (10) reduce to a single one 



F q [h] 



dh(q;z) 



z)-c(q;s)} (14) 



where F q [h] 



c(q;s) - h(q;s) + p[h(q; s) Ac{q) + 
c(q;s)Ah(q) + c(q;0)h{q; s) - {h{q; s) - h{q;0))c(q; s)}. 
This equation is a particular MCT-like equation whose 
memory kernel is represented by the direct correlation 
function extracted from 



hi[h(x; s) + 1] = h{x: s) — c(x; s) 



(15) 



The phenomenology displayed by the above dynami- 
cal equation is similar to the one of conventional MCT 
and its solution can be identified with the scaling func- 
tion determining the correlation function within the 
temperature-time superposition principle. The equation 
is invariant under time rescaling t — > At for positive A, 
and the condition for existence of a decaying solution is 



that the operator 



dF q [h] 



has a zero mode K.o(q) which 



dh(k) 

has been computed in [7, 8]. By analyzing the structure 
of the equation one can derive an analytic expression for 
the MCT exponent parameter 



A = 



fd D xK%(x)/(l + h(,x))* 



pAc(g)]3 



(16) 



In [7, 8] this exponent parameter has been computed us- 
ing a completely different static approach. Here we want 
to stress that our result can be obtained using a full dy- 
namical description. The a relaxation time r a can also 
be extracted from the equation, through standard MCT 
analysis and one find r Q ~ \T — Ta\ ' with 7 = w- + 55. 

Aging solution. Below the dynamical glass transition, 
the dynamical equations (9) and (10) describe an ag- 
ing regime, under the assumption of a high temperature 
initial condition and loss of memory of the initial con- 
dition h(q;s,0) = c(g;s,0) = 0. The resulting equa- 
tions can be solved using the modified FDT relation 
Rh{q;t,s) = X/3dh(q;t,s)/ds [14]. As usual, the com- 
bination T/X is interpreted as an effective temperature 
[9], and coherently X is independent of q. As it has been 
shown in [14], the FDT ratio X E [0, 1] is fixed by the re- 
quirement that the equation for the response function in 
the limit u — ^ s admits a non vanishing solution. Again 
this is only possible in presence of a zero mode ICq (q) of 
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FIG. 1: The value of X as a function of ho in the schematic 
model for p — 0.8. We neglect all space dependence in the 
equations and seek for a solution of the equations as a function 
of ho. Higher values of ho represent stronger attraction. For 
ho = 1.403 there is a dynamical glass transition. The FD- 
ratio is 1 at the transition and tends monotonically to for 
ho — > 00. Inset the value of Ah as a function of ho- 



the kernel operator: 
6(q, k) = (2pAc{q) - p 2 Ac 2 (q)) (2n) D 5(q - k)- 



I 



1 + h 



(q-k) 



(17) 



Examples of quantities that can be computed with our 
approach are presented in figure 1. Instead of solving the 
full dynamical equations for three dimensional liquids wc 
have simplified them in correspondence of the first peak 
of the static structure function Iiq = limt^.ooh(t,t, q pea k) 
taken as input and neglected any momentum and space 
dependence. The resulting values of X and h are plotted 
as functions of ho for p = 0.8. 

A different closure: conventional MCT equations. The 
crucial point to obtain the dynamical equations (9-14) 
has been that, modulo a different interpretation of the 
replica index, there is formal coincidence of replica pseu- 
dodynamic HNC equations with the ordinary equations 
of replicated liquids. This coincidence, which can be 
traced in the replica symmetry of the pseudodynamic 
effective Hamiltonian, should be respected in any con- 
sistent approximation scheme describing the correlation 
functions. It has been recently observed in [5] that the 
MCT equation specifying the non-ergodicity parameter 
can be derived from replicated Ornstein-Zernike equa- 
tion, through a defined closure of Yvonne-Born-Green 
(YBG) hierarchy leading to the following approximation 
for the non diagonal elements of the replica directed cor- 
relation function c a p(q) for a f3 



ap(k) = dq V(k,q)h a/3 (q)h al3 (k - q) 



(18) 
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where the V(k, q) is the Mode Coupling vertex given in 
terms of the static direct correlation function cq (q) as 



V(k,q) = 



1 



167T 3 /fc 2 



[k.(q CQ ( (Z ) + (k-q)c (fc-g))] 2 ,(19) 



which is independent of the replica indexes. Consis- 
tently with our previous comment, we can then inter- 
pret (18) in pscudodynamics with a — > (a = l,t), 
j3 (a = 1, s = 0) with t > 0. Assuming TTI, we 
get c(k,t) = J dq V(k,q)h(q,t)h(k — q,t). Inserting in 
(14) one easily see that exactly the critical Gotze mode 
coupling equation is obtained. We notice that in [17] 
MCT has been generalized to describe aging below T<j. 
We leave for future work the analysis of possible gener- 
alizations Szamcl closure to this case. 
Conclusions. Summarizing, we have shown that within 
liquid theory approximations, under marginal conditions, 
Boltzmann pscudynamics gives rise to a dynamical pic- 
ture strictly following the MCT phenomenology. We have 
been able to derive dynamical MCT-like equations in 
a coherent scheme that includes both equilibrium and 
aging dynamics in the context of HNC approximation. 
These equations allow to describe equilibrium dynam- 
ics on approaching and aging dynamics below 
We have computed the zero mode and the A param- 
eter from these equations. In equilibrium both these 
quantities coincide with the ones computed with com- 
pletely static techniques in [7]. In aging, the descrip- 
tion is coherent with the p-spin picture, but goes beyond 
this schematic model. In this case we have derived an 
expression for the FDT ratio. We also found that the 
conventional Gotze MCT equation can be obtained from 
our replica-dynamic Ornstein-Zernikc equations general- 
izing to pseudo-dynamics a closure scheme of YBG hi- 
erarchy proposed by Szamel. Our analysis demonstrates 
that critical MC and MC-likc dynamical behavior in equi- 
librium and aging is equivalent to quasi-equilibrium sam- 
pling of configuration space. We expect a similar picture 
to hold in other approximation schemes. In particular it 
will be interesting to study the case of small cage approx- 
imation that has been found to give exact result in the 
limit of high dimensionality [18]. These are mean- field 
like results. The usual caveats that apply to mean- field 
theories are valid here. Activated processes are neglected 
and the dynamical transition is spurious. To our view, 
quasi-equilibrium exploration of configuration space is a 
fundamental property of glassy dynamics that should also 
hold when mean- field theory fails. Theoretical and exper- 
imental efforts will be needed to test this conjecture. 
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